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PREFACE 


In  March  1987,  the  Norfolk  District  Corps  of  Engineers  requested  that  the 
US  Army  Engineer  Waterways  Experiment  Station  (WES)  conduct  an  investigation  to 
assess  general  changes  in  circulation,  currents,  and  sedimentation  associated 
with  six  proposed  alternative  expansion  geometries  of  Craney  Island,  the  con¬ 
fined  dredged  material  disposal  site  located  in  the  lower  James  River. 

The  study  was  conducted  by  personnel  of  the  Hydraulics  Laboratory,  WES, 
under  the  general  direction  of  Messrs.  F.  A.  Herrmann,  Jr.,  Chief  of  the 
Hydraulics  Laboratory;  R.  A.  Sager,  Assistant  Chief  of  the  Hydraulics  Labora¬ 
tory;  W.  H.  McAnally,  Jr.,  Chief  of  the  Estuaries  Division;  and  W.  D.  Martin, 
Chief  of  the  Estuarine  Engineering  Branch.  The  project  was  conducted  by 
Messrs.  S.  B.  Heltzel  and  M.  A.  Granat,  Estuarine  Engineering  Branch.  This 
report  was  prepared  by  Messrs.  Heltzel  and  Granat  and  edited  by  Mrs.  Marsha  C. 
Gay  of  the  Information  Technology  Laboratory,  WES. 

The  valuable  technical  and  nontechnical  contributions  of  Mr.  J.  R. 
Melchor,  Norfolk  District  Corps  of  Engineers,  are  gratefully  acknowledged. 

COL  Dwayne  G.  Lee,  CE,  is  the  Commander  and  Director  of  WES. 

Dr.  Robert  W.  Whalin  is  the  Technical  Director. 
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CONVERSION  FACTORS,  NON-SI  TO  ST  (METRIC) 
UNITS  OF  MEASUREMENTS 


Non-SI  units  of  measurement  used  in  this  report  can  be  converted  to  SI  (metric) 
units  as  follows: 


Multiply 

acres 

cubic  feet 
cubic  yards 
feet 


_ gZ  , 

4,046.873 

0.02831685 

0.7645549 

0.3048 


_ To  Obtain 

square  metres 
cubic  metres 
cubic  metres 
metres 
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LOWER  JAMES  RIVER  CIRCULATION  STUDY,  VIRGINIA 
EVALUATION  OF  CRANEY  ISLAND  ENLARGEMENT  ALTERNATIVES 


PART  I:  INTRODUCTION 

Background 

1.  The  existing  Craney  Island  disposal  area  is  a  2,500-acre*  confined 
disposal  facility  located  near  Norfolk,  Virginia  (Figure  1).  Design  plans  for 
the  facility  were  developed  in  the  early  1940's,  and  construction  was  under¬ 
taken  between  1954  and  1957.  Several  retaining  dike  elevation  increases  and 
the  implementation  of  a  management  plan  developed  specifically  for  Craney 
Island  (Palermo,  Shields,  and  Hayes  1981)  have  greatly  extended  the  useful 
life  of  the  facility  beyond  initial  expectations.  At  present,  approximately 

5  million  cubic  yards  of  predominately  fine-grained  maintenance  dredged  mate¬ 
rial  from  the  channels  and  ports  in  the  Hampton  Roads  area  are  annually  dis¬ 
posed  of  in  the  facility. 

2.  Additional  new  work  dredged  material  from  the  Norfolk  Harbor  and 
channels  deepening  and  improvement  project  will  also  be  placed  into  Craney 
Island.  This  new  work  dredged  material  will  greatly  reduce  the  storage  capac¬ 
ity  aud  active  life  of  Craney  Island.  An  expansion  of  the  existing  site  or 
the  development  of  an  alternate  site  will  be  necessary  to  contain  future  new 
work  and  maintenance  dredged  material  once  the  existing  disposal  site  has 
reached  its  capacity.  Six  alternative  expansion  configurations  for  Craney 
Island  are  presently  being  considered  by  the  Norfolk  District  Corps  of 
Engineers . 

3.  An  assessment  of  the  impact  of  each  alternative  on  circulation  and 
sedimentation  in  the  lower  James  River  was  undertaken  as  a  preliminary  plan¬ 
ning  level  task.  Other  concurrent  tasks  included  analyses  to  determine  the 
storage  capacity  and  active  dredged  material  disposal  life  of  each  of  the 
proposed  alternatives  and  development  of  preliminary  guidelines  for  management 


*  A  table  of  factors  for  converting  non-SI  to  SI  (metric)  units  of  measure¬ 
ment  is  presented  on  page  3. 
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of  the  new  areas*  and  an  evaluation  of  the  required  retaining  dike  geotech¬ 
nical  engineering  characteristics  (stability,  constructibility ,  and  cost)  for 
each  alternative  (Spigolon  and  Fowler  1987). 

Purpose 

4.  The  primary  objective  of  the  present  study  was  to  use  available 
numerical  models  to  assess  general  changes  in  circulation,  currents,  and  sed¬ 
imentation  associated  with  each  of  the  six  proposed  alternative  expansion 
geometries  of  Craney  Island.  An  additional  objective  of  the  study  was  to 
assess  the  effects  of  each  of  the  six  alternative  geometries  on  the  reported 
estuarine  circulation  cell  (flow  convergence)  off  Hampton  Flats  and  Newport 
News  Point. 


Scope 

5.  The  numerical  modeling  portion  of  this  study  was  designed  to  address 
relative  alternative-induced  changes  in  overall  hydrodynamics  and  to  assess 
relative  sedimentation  changes  in  four  specific  critical  zones  of  interest. 

The  primary  focus  was  circulation  and  sedimentation  in  relatively  low  velocity 
areas,  so  the  cohesive  version  of  a  numerical  sediment  code  was  used  since  it 
more  directly  reflected  circulation  and  suspended  sediment  transport.  Fig¬ 
ure  2  illustrates  the  James  River  study  area  with  the  four  zones  of  interest 
highlighted:  Area  A,  the  entrance  to  Willoughby  Bay;  Area  B,  the  Hampton 
Flats;  Area  C,  the  entrance  to  the  Nansemond  River;  and  Area  D,  Burwell  Bay. 

6.  The  circulation  cell  off  Newport  News  Point  and  Hampton  Flats  was 
addressed  using  results  from  the  numerical  model  study,  comparison  with  sur¬ 
face  current  pattern  mosaics  from  earlier  physical  model  investigations** 

of  the  James  River,  review  of  a  recent  investigation  of  the  area  (Byrne  et  al. 
1987)  conducted  by  the  Virginia  Institute  of  Marine  Science  (VIMS) ,  and  a 


*  G.  F.  Goforth.  1986.  "Disposal  Life  Evaluation  of  Alternative  Expansion 
Configurations  for  Craney  Island  Disposal  Facility,"  Draft  Report,  US  Army 
Engineer  Waterways  Experiment  Station,  Vicksburg,  Miss. 

**  N.  J.  Brogdnn ,  Jr.,  and  W.  H.  Bobb.  1967.  "Effects  of  Proposed  Water¬ 
front  Developments  at  Newport  News  Point  on  Tides,  Currents,  Salinities, 
and  Shoaling."  Draft  Renort,  US  Army  Engineer  Waterways  Experiment  Station, 
Vicksburg,  Miss. 
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review*  of  a  limited  field  data  set  collected  by  the  US  Army  Engineer  Water¬ 
ways  Experiment  Station  (WES)  during  a  relatively  low  freshwater  discharge 
period  during  July  1986. 
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*  S.  B.  Heltzel.  1986  (2  Sep).  "Memorandum  for  Record:  Interim  Report  - 
Norfolk  Harbor  Long-Term  Disposal  Study,"  US  Army  Engineer  Waterways 
Experiment  Station,  Vicksburg,  Miss. 
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PART  II:  NUMERICAL  MODELING  APPROACH 


!•- 


The  Numerical  Models 


7.  The  Corps  numerical  modeling  system,  Open-Channel  Flow  and  Sedimen¬ 
tation,  TABS-2  (Thomas  and  McAnally  1985),  was  used  in  this  investigation. 

The  two  primary  finite  element  numerical  model  codes  used  were  A  Two- 
Dimensional  Model  for  Free  Surface  Flows  (RMA-2V)  and  Sediment  Transport  in 
Unsteady  Two-Dimensional  Flows,  Horizontal  Plane  (STUDH) .  Both  codes  employ 
the  finite  element  method  to  solve  the  depth-integrated  governing  equations. 
Appendix  A  contains  general  information  on  the  finite  element  method.  A  brief 
description  of  RMA-2V  and  STUDH  appears  in  Appendices  B  and  C,  respectively. 

Lower  James  River  Computational  Meshes 

8.  The  computational  mesh  used  during  this  investigation  was  a  modified 
version  of  the  mesh  developed  for  the  1-664  bridge  tunnel  crossing  study 
(Heltzel ,  in  preparation).  Modifications  to  this  mesh  included  (a)  revising 
the  schematization  around  Craney  Island  by  adding  additional  elements  to 
incorporate  the  six  Craney  Island  enlargement  alternatives,  (b)  extending  the 
mesh  into  lower  Chesapeake  Bay  to  improve  modeling  of  hydrodynamic  and  sedi¬ 
mentation  processes  in  the  Willoughby  Bay  area,  and  (c)  effectively  doubling 
the  mesh  resolution  in  the  critical  zones  of  Burwell  Bay,  the  Nansemond  River 
entrance,  and  Newport  News  Point/Hampton  Flats  area.  All  of  the  conditions 
tested  included  the  completed  1-664  bridge  tunnel  crossing  and  the  55-ft  New¬ 
port  News  and  Norfolk  Harbor  channels. 

9.  The  basic  revised  mesh,  presented  in  Figure  3,  contains  2,326  nodes 
and  806  elements.  This  mesh  incorporates  all  six  alternative  geometries  to 
eliminate  the  possibility  of  required  mesh  resolution  refinement  between  con¬ 
ditions  as  a  possible  cause  for  anomalous  plan  variations  in  hydrodynamic  or 
sedimentation  results.  An  enlarged  view  of  the  Craney  Island  area  illustrat¬ 
ing  the  schematization  used  for  each  of  the  proposed  enlargement  alternatives 
is  provided  in  Figure  4.  The  highlighted  areas  in  each  schematization  indi¬ 
cate  the  elements  that  were  deleted  from  the  computational  space  for  each  of 
the  conditions  considered.  During  testing,  the  boundary  of  each  mesh  con¬ 
formed  to  the  new  geometry  for  each  of  the  respective  alternatives.  Figure  5 
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Plan  A 


Plan  B 


Plan  C 


a.  Plan  D 


Plan  E 


Plan  F 


Figure  5.  Dimensions  of  Craney 
Island  enlargement  alternatives 


provides  the  actual  plan  dimensions  for  each  of  the  analyzed  Craney  Island 
extensions . 


Testing  Conditions 


10.  Physical  model  data  collected  in  the  Chesapeake  Bay  hydraulic  model 
during  the  Norfolk  Harbor  and  channels  deepening  study  (Richards  and  Morton 
1983)  were  used  to  develop  boundary  conditions  for  the  numerical  hydrodynamic 
model.  Physical  model  water-surface  elevation  data  were  used  to  generate 
boundary  forcing  functions  at  the  lower  Chesapeake  Bay  numerical  model  bound¬ 
ary,  and  depth-averaged  physical  model  velocity  data  were  used  at  the  upper 
James  and  Elizabeth  River  numerical  model  boundaries. 

11.  The  same  hydrodynamic  coefficients  and  modeling  procedures  devel¬ 
oped  during  the  1-664  bridge  tunnel  crossing  study  (Heltzel,  in  preparation) 
and  the  numerical  Norfolk  Harbor  and  channels  deepening  study  (Berger  et  al. 
1985)  were  used  during  this  investigation.  Water-surface  elevation  and  veloc¬ 
ity  data  from  several  interior  locations  within  the  revised  mesh  were  compared 
to  data  from  the  1-664  investigation  to  ensure  that  the  revised  schematization 
and  boundary  conditions  did  not  alter  hydrodynamic  characteristics.  Cohesive 
sediment  coefficients  and  modeling  procedures  developed  during  the  1-664  and 
Norfolk  investigations  were  similarly  employed  during  the  Craney  Island  study. 
These  parameters  provided  results  that  compared  well  with  the  limited  field 
sedimentation  information  and  the  suspended  sediment  concentrations  reported 
for  the  areas  of  interest. 

12.  A  series  of  base  (revised  mesh  without  any  Craney  Island  expansion) 
numerical  model  sensitivity  studies  were  undertaken  prior  to  final  base  and 
plan  testing  to  determine  the  boundary  condition  that  produced  the  maximum 
sedimentation  in  the  areas  of  interest.  As  indicated  in  the  numerical  Norfolk 
Harbor  and  channels  deepening  study  (Berger  et  al.  1985),  and  as  confirmed  by 
these  sensitivity  studies,  the  mean  range  tide  (2.5  ft  at  Old  Point  Comfort) 
and  long-term  average  James  River  freshwater  discharge  condition  (8,900  cfs 
combined  total  James  River  tributary  freshwater  inflow)  generally  resulted  in 
maximum  sedimentation  rates.  The  mean  range  tide  and  long-term  average  fresh¬ 
water  discharge  conditions  were  used  as  the  forcing  functions  for  final  test¬ 
ing  of  the  base  and  plan  configurations.  It  should  be  noted  that  the  ocean 
salinity  was  maintained  at  32.5  ppt  during  the  physical  model  study. 
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PART  III:  ESTUARINE  CIRCULATION  AND  FLOW  CONVERGENCE: 
HAMPTON  FLATS  AND  NEWPORT  NEWS  POINT 


13.  This  section  summarizes  information  presently  available  regarding 
the  estuarine  circulation  and  flow  convergence  observed  off  Newport  News  Point 
and  the  Hampton  Flats.  This  summary  is  based  on  previous  studies*  conducted 
on  the  James  River  physical  model  at  WES,  Vicksburg,  Mississippi;  a  recent 
detailed  investigation  conducted  by  VIMS  (Byrne  et  al.  1987)  evaluating  poten¬ 
tial  impacts  associated  with  the  development  of  an  island  (New  Port  Island)  on 
Hampton  Flats;  and  a  limited  supplemental  field  data  collection  effort  con¬ 
ducted  by  WES  during  the  period  22-24  July  1986.** 

James  River  Physical  Model  Investigations 

14.  Several  physical  model  dye  studies  in  the  late  1960's  (as  sum¬ 
marized  in  Byrne  et  al.  1987)  investigated  James  River  circulation  patterns  to 
determine  optimum  oyster  larvae  release  locations  in  an  attempt  to  reestablish 
the  lower  James  River  oyster  beds.  These  early  investigations  demonstrated 
the  existence  and  importance  of  the  estuarine  net  nontidal  flow  character¬ 
istics  of  the  lower  James  River. 

15.  Surface  current  pattern  mosaics  taken  during  several  of  the  physi¬ 
cal  model  studies  provide  excellent  visual  documentation  of  the  existing  sur¬ 
face  circulation  patterns.  Figure  6  illustrates  a  set  of  photographs  taken 
over  a  4-hour  period  (2.4  min  in  the  model)  during  a  study  sponsored  by  the 
City  of  Newport  News  to  specifically  investigate  effects  of  proposed  water¬ 
front  developments  at  Newport  News  Point  on  tides,  currents,  salinities,  and 
shoaling.  The  study  was  conducted  during  the  period  November  1966  to  March 
1967,  and  the  results  were  documented  by  Brogdon  and  Bobb  (1967).* 

16.  Each  photograph  illustrates  the  trajectory  of  surface  confetti 
during  a  simulated  5-min  period  (3  sec  in  the  model) .  A  strobe  light  acti¬ 
vated  at  the  end  of  each  exposure  dotted  the  confetti  streaks,  identifying  the 
flow  direction.  The  indicated  times  are  referenced  to  the  specific  hour  after 
the  moon's  transit  over  the  entrance  to  Chesapeake  Bay  (the  76th  meridian). 

*  Brogdon  and  Bobb,  op.  cit. 

**  He  1 tze 1 ,  op .  cit. 
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The  illustrated  base  condition  included  the  proposed  35-ft  channel  from  New¬ 
port  News  to  Richmond  and  the  proposed  enlargement  of  Newport  News  Channel  and 
anchorages.  The  boundary  conditions  included  a  2.5-ft  mean  tide  range  condi¬ 
tion  at  Hampton  Roads,  a  Chesapeake  Bay  sump  salinity  of  24.2  ppt,  and  the 
long-term  average  freshwater  inflow  of  7,500  cfs  for  the  James  at  Richmond, 
1,000  cfs  for  the  Appomattox,  300  cfs  for  the  Chickahominy ,  and  700  cfs  for 
the  Nansemond  which  included  the  discharges  of  the  Warwick,  Pagan,  Chuckatuck, 
and  Elizabeth  rivers. 

17.  As  indicated  in  Figure  6,  an  early  onset  of  the  flood  current  over 
the  Hampton  Flats  area  began  at  hour  3  while  the  flow  upstream  of  Newport  News 
Point  and  south  of  Hampton  Flats  was  still  at  strength  of  ebb.  As  demon¬ 
strated  at  hour  4,  flood  flow  over  Hampton  Flats  was  fully  developed  while 
flows  west  of  Newport  News  Point  and  in  Newport  News  Channel  were  still 
ebbing.  The  flood  current  was  not  fully  developed  in  Newport  News  Channel  and 
we«t  of  Newport  News  Point  until  hour  6.  Similar  circulation  patterns  were 
demonstrated  in  the  other  physical  model  studies. 

New  Port  Island  Investigation 

18.  The  New  Port  Island  study  recently  completed  by  VIMS  (Byrne  et  al. 
1987)  was  an  extensive  multitasked  investigation  conducted  to  evaluate  poten¬ 
tial  lower  James  River  marine  resources  impacts  associated  with  the  develop¬ 
ment  of  New  Port  Island,  a  recreational/port  facility,  to  be  located  on 
Hampton  Flats.  A  major  emphasis  of  this  study  was  the  estuarine  circulation 
off  Newport  News  Point,  its  effects  on  oyster  larvae  transport,  and  the  poten¬ 
tial  impacts  associated  with  various  alternative  configurations  for  the  pro¬ 
posed  island. 

19.  The  VIMS  study  concluded  that  the  flow  convergence  off  Newport  News 
Point,  described  as  a  frontal  system,  was  a  persistent  phenomenon  expected  to 


occur  during  times  of  flood  current  under  normal  tidal  range,  freshwater 
inflow,  and  meteorological  conditions.  The  flow  convergence  was  said  to  be 
the  result  of  geometry-induced  local  phase  (time  of  arrival)  differences 
between  initial  flood  currents  on  Hampton  Flats  and  Newport  News  Channel  and 
ebbing  flows  west  of  Newport  News  Point.  The  local  salinity  distribution  was 
said  to  be  another  essential  factor  to  the  formation  of  the  front.  During 
development  of  the  flow  convergence,  the  currents  on  Hampton  Flats,  and  later 
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in  the  Newport  News  Channel,  began  to  flood  while  the  currents  west  of  Newport 
News  Point  were  still  ebbing.  The  study  also  concluded  that  the  depth  tran¬ 
sition  at  Newport  News  Point  enhanced  and  stabilized  the  location  of  the  con¬ 
vergence  zone.  Dye  studies  demonstrated  that  the  surface  floodwaters 
downstream  of  the  convergence  zone  were  transported  to  depths  of  13  to  16  ft 
upstream  of  the  convergence  zone. 

WES  Supplemental  Field  Survey 

20.  WES  conducted  a  limited  field  survey*  in  the  lower  James  River  on 
22-24  July  1986  to  acquire  additional  supplemental  field  data  during  a  pre¬ 
liminary  stage  in  the  development  of  a  three-dimensional  James  River  modeling 
effort.  Four  boats  were  used  to  collect  vertical  velocity  and  salinity  data 
at  13  stations  during  a  complete  13-hr  tidal  cycle  on  22  and  24  July.  Data 
were  collected  at  five  locations  in  the  vertical  and  at  approximately  half- 
hour  intervals  where  station  spacing  permitted. 

21.  Environmental  conditions  were  not  representative  of  normal  condi¬ 
tions  in  the  lower  James  River.  The  freshwater  discharge  of  approximately 
1,000  cfs  was  well  below  the  long-term  average  condition.  In  addition,  strong 
winds  were  blowing  and  gusting  from  the  east  and  southeast  during  the  survey 
effort.  The  boat  located  at  Newport  News  Point  was  unable  to  complete  the 
entire  tidal  cycle  survey  due  to  high  wave  conditions;  data  were  collected 
only  during  the  later  stages  of  ebb  and  the  early  stages  of  flood,  the  pre¬ 
dicted  period  of  time  for  the  convergence  phenomenon.  As  confirmed  by  the 
vertical  salinity  data,  these  conditions  resulted  in  a  reduced  stratification. 

22.  Instantaneous  middepth  velocities  were  analyzed  for  the  13  stations 
during  this  period  of  time.  These  data  demonstrated  the  early  flood  on 
Hampton  Flats  relative  to  the  Newport  News  Channel  station.  Tidally  averaged 
bottom  velocities  generally  illustrated  net  upstream  movement  onto  the  Hampton 
Flats  and  net  downstream  movement  in  the  Newport  News  Channel.  The  velocities 
averaged  over  depth  and  then  over  the  tidal  cycle  illustrated  the  same  general 
trends  with  some  magnitude  changes.  The  existence  of  this  circulation  and  the 
convergence  phenomenon,  despite  the  unfavorable  conditions  for  its  formation 
'low  freshwater  discharge  and  high  winds  and  waves),  clearly  demonstrates  its 


/ 1*  1  ,  op  .  c  i  t  . 
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PART  IV:  MODELING  RESULTS  AND  DISCUSSION 


Hydrodynamic  Impacts 

23.  Localized  and  subtle  base  to  plan  hydrodynamic  differences  were 
indicated.  Hydrodynamic  vector  plots  for  maximum  ebb*  (hr  14.00)  and  maximum 
flood  (hr  20.00)  conditions  in  the  lower  James  River  for  each  of  the  plan  con¬ 
ditions  are  printed  in  red  in  Plates  1-12.  The  corresponding  base  condition 
is  printed  in  black  on  each  plate.  Except  for  the  Craney  Island  extension 
area,  where  only  base  vectors  occur,  the  plan  velocity  was  identical  to  base 
conditions  (i.e.,  when  red  plan  vectors  exactly  overlay  black  base  vectors, 
only  the  black  vector  is  visible)  .  Closeup  hydrodynamic  vector  plots  for  the 
Newport  News/Hampton  Flats  area  for  the  period  around  slack  before  flood 

(hr  15.00,  16.00,  and  17.00)  for  each  of  the  plan  conditions  are  printed  in 
red  in  Plates  13-18.  These  plots  illustrate  the  Hampton  Flats  circulation 
cell  and  the  convergence  zone  for  each  plan  condition.  Again,  the  correspond¬ 
ing  base  circulation  is  printed  in  black  on  each  plot.  Each  of  the  vector 
plots  is  based  on  a  regular  grid  pattern  that  uses  the  finite  element  shape 
functions  and  the  calculated  nodal  velocity  vectors. 

24.  Subtle  localized  variations,  generally  within  16,000  ft  adjacent  to 
and  north  and  northwest  of  Craney  Island,  are  indicated  on  the  vector  plots 
comparing  the  plan  conditions  to  the  base  condition.  Plans  A-C  (Plates  1-6), 
which  all  involve  northward  Craney  Island  extensions,  appear  to  illustrate  the 
largest  hydrodynamic  impacts.  Plans  A  and  B,  which  also  involve  westward 
extensions,  illustrate  the  largest  changes.  Some  of  the  indicated  variation 
may  be  the  result  of  small  plan  to  base  phase  shifts  or  simply  numerical 
noise. 

25.  The  vector  plots  provide  an  excellent  visual  presentation  of  the 
circulation  for  each  condition;  however,  quantifying  variations  between  the 
various  plans  is  difficult  from  these  plots.  Figure  7  illustrates  the  loca¬ 
tions  of  13  nodes  that  were  examined  in  detail  to  summarize  the  actual 

*  As  in  the  physical  model  investigations,  the  indicated  times  are  referenced 
to  the  moon's  transit  over  the  entrance  to  Chesapeake  Bay  (the  76th  merid¬ 
ian).  A  repetitive  !2.42-hr  tidal  cycle  was  used  in  the  numerical  model 
runs;  i.e.,  hr  13.00  (the  first  time-step  in  cycle  2)  corresponds  to  hr  0.58 
of  the  first  tidal  cycle. 


differences  in  magnitude  between  the  velocities  of  the  base  condition  and  each 
of  the  plan  conditions.  In  the  order  listed,  node  18  is  located  in  the  center 
of  the  Willoughbv  Bay  zone  of  interest  (Area  A);  nodes  90,  146,  and  198  are 
located  in  Hampton  Flats  (Area  B) ;  nodes  598,  677,  and  2118  are  located  in  the 
Nansemond  River  zone  of  interest  (Area  C) ;  nodes  1218,  1258,  and  1298  are 
located  in  the  Burwell  Bay  zone  of  interest  (Area  D) ;  and  nodes  1994,  1980, 
and  1972  are  located  in  the  thalweg  of  the  Newport  News  Channel. 

26.  Table  1  provides  the  maximum  ebb  and  flood  magnitudes  for  base  con¬ 
ditions  at  each  node  location  and  the  magnitude  differences  for  each  of  the 
plan  conditions.  As  indicated,  with  the  exception  of  Newport  News  Channel, 
all  maximum  velocity  differences  (plan  magnitude  minus  base  magnitude)  at  the 
examined  critical  areas  of  interest  were  generally  within  field  and  model  con¬ 
fidence  limits  and  never  greater  than  0.06  fps  from  base  conditions.  Notice¬ 
able  plan  to  base  differences  (velocity  differences  greater  than  0.10  fps) 
were  indicated  for  the  Newport  News  Channel  nodes.  Channel  plan  velocities 
always  exceeded  base  velocities  with  maximum  ebb  velocity  differences  greater 
than  maximum  flood  velocity  differences.  As  expected,  the  plans  with  north¬ 
ward  extensions  resulted  in  the  largest  increases.  The  greatest  changes,  less 
than  0.35  fps  on  ebb  and  0.25  fps  on  flood  (Table  1),  were  indicated  for  plans 
A  and  B,  the  largest  expansion  alternatives  which  also  involved  westward 
expansions. 

27.  Plates  19-28  illustrate  selected  node  time-history  plots  of  water- 
surface  elevations  and  ebb  and  flood  velocity  magnitudes.  A  node  from  the 
center  of  each  of  the  critical  areas  and  one  from  the  Newport  News  Channel  was 
selected  from  each  plan  condition  for  comparison  to  the  base  condition.  Sim¬ 
ilar  plots  for  other  node  locations  are  available,  but  are  not  included  in 
this  report.  The  first  few  hours  of  each  model  run  should  not  be  used  for 
comparison  or  analysis  since  they  may  be  affected  by  model  spin-up  conditions. 

28.  The  time-history  plots  provide  an  excellent  means  of  illustrating 
actual  phase  (time  of  arrival)  and  magnitude  differences  at  specific  node 
locations  over  the  tidal  cycle.  As  previously  addressed,  with  the  exception 
of  the  Newport  News  Channel,  only  subtle  base  to  plan  hydrodynamic  variations 
were  identified.  The  subtle  phase  shifts  illustrated  may  be  partially  respon¬ 
sible  for  the  some  of  the  variations  indicated  on  the  vector  plots. 

29.  Results  from  the  depth-integrated  numerical  model  demonstrate  that 
the  formation  of  the  circulation  cell  off  Newport  News  and  Hampton  Flats  will 


continue  under  all  of  the  alternative  plans  examined.  This  circulation  cell 
is  the  result  of  the  early  onset  of  the  flood  currents  over  Hampton  Flats 
associated  with  bathymetric  shielding  of  the  ebb  currents  over  the  flats  by 
Newport  News  Point.  This  large-scale  geometry-induced  phase  difference 
existed  during  all  of  the  plans  tested. 

30.  The  differences  indicated  in  the  time-history  and  vector  plots  are 
subtle  and  indicate  the  impacts  to  be  localized,  within  16,000  ft  of  Craney 
Island  (see  paragraph  24),  with  no  impacts  identified  to  the  general  estuarine 
circulation  outside  of  this  region.  The  bathymetric  shielding,  combined  with 


the  depth  transition  as  described  in  the  VIMS  study  (Byrne  et  al .  1987)  and 
summarized  in  paragraph  19,  has  the  effect  of  stabilizing  the  location  of  the 
convergence  zone. 

31.  Potential  impacts  to  the  reported  frontal  system,  a  three- 
dimensional  phenomenon,  should  not  be  quantified  with  the  information  from  the 
present  investigation.  The  following  generalizations  can  be  made,  however, 
based  upon  the  available  information  and  present  estuarine  experience.  Hydro- 
dynamic  impacts  to  the  general  three-dimensional  circulation  will  be  small,  if 
any,  and  will  be  extremely  difficult  to  measure.  Alternatives  with  northward 
Craney  Island  expansions  will  have  the  greatest  potential  to  affect  the 
three-dimensional  circulation  processes. 

Sedimentation  Impacts 

32.  Sedimentation  rates  are  generally  sensitive  to  subtle  variations 
in  circulation  characteristics.  Small  changes  in  hydrodynamic  processes  are 
usually  amplified  in  the  sedimentation  responses.  Sedimentation  comparisons 
therefore  provide  an  excellent  means  of  quantitatively  assessing  overall  base 
to  plan  impacts. 

33.  For  reporting  purposes,  the  predicted  shoaling  volume  for  each 
element  in  the  zone  of  interest  was  combined  to  estimate  yearly  shoaling  rates 
for  each  zone.  Table  2  provides  a  summary  of  the  plan-predicted  sedimentation 
divided  by  the  base-predicted  sedimentation,  the  plan  to  base  shoaling  index 
value,  for  the  four  critical  areas  of  interest  for  each  of  the  six  alternative 
plans.  Based  on  available  field  data,  shoaling  rates  in  each  of  these  areas 
are  generally  low  (less  than  0.5  ft/year),  and  the  indicated  shoaling  index 
values  are  generally  well  within  normal  hydrographic  survey  detection  limits. 
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34.  Of  the  four  zones  examined,  the  Nansemond  River  entrance,  Area  C, 
is  the  only  zone  considered  to  demonstrate  any  distinct  change  in  base  to  plan 
sedimentation.  It  should  be  stressed  that  even  these  variations  are  rather 
subtle,  especially  considering  the  generally  low  sedimentation  rates. 

Plans  A,  B,  and  C,  which  all  involve  northward  extensions  of  Craney  Island, 
resulted  in  reduced  sedimentation  in  the  Nansemond  area  for  the  alternative 
conditions.  The  largest  sedimentation  change,  about  a  9  percent  reduction  in 
shoaling,  was  indicated  for  plan  C,  the  alternative  that  involved  expansion 
only  to  the  north. 

35.  Plans  D,  E,  and  F,  which  involve  only  westward  extensions  of  Craney 
Island,  resulted  in  increased  sedimentation  in  the  Nansemond  area.  Plan  F 
resulted  in  the  largest  increase,  about  a  7  percent  increase  from  the  base 
shoaling  rate.  The  second  largest  increase,  about  4  percent,  was  indicated 
for  plan  D.  It  is  interesting  to  note  that  none  of  alternatives  A,  D,  nor  F 
extends  to  the  mainland  (Figure  4) ,  allowing  additional  sluggish  circulation 
between  the  disposal  site  and  the  mainland.  Plan  E  extends  to  the  mainland. 
Although  this  alternative  encompasses  a  larger  area  than  either  plans  F  or  D, 
it  appears  to  have  a  reduced  impact  on  sedimentation.  An  additional  consid¬ 
eration  associated  with  plans  A,  D  and  F  is  a  potential  impact  on  water 
quality  associated  with  reduced  circulation  between  the  mainland  and  the  dis¬ 
posal  site  (Plates  1  and  2,  7  and  8,  and  11  and  12). 
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PART  V:  CONCLUSIONS 


36.  No  plan  to  base  velocity  differences  greater  than  ±0.06  fps  were 
identified  at  any  of  the  four  critical  areas  of  interest.  With  the  exception 
of  the  Newport  News  Channel,  only  localized  and  subtle  base  to  plan  velocity 
differences  were  indicated.  Plan  to  base  magnitude  increases  greater  than 
±0.10  fps  were  identified  in  the  Newport  News  Channel.  The  greatest 
increases,  less  than  0.35  fps  on  ebb  and  0.25  fps  on  flood,  were  indicated  for 
plans  A  and  B,  which  involved  both  northward  and  westward  Craney  Island 
extensions . 

37.  Subtle  localized  circulation  variations  were  identified  in  base  to 
plan  comparison  vector  plots.  When  these  variations  occurred,  they  were  gen¬ 
erally  within  16,000  ft  north  and  northwest  of  Craney  Island.  In  some 
instances,  these  variations  were  associated  with  slight  hydrodynamic  phase 
shifts. 

38.  Plan  to  base  shoaling  index  values  (plan-predicted  sedimentation 
divided  by  base-predicted  sedimentation)  for  all  conditions  were  well  within 
hydrographic  survey  detection  limits.  The  Nansemond  River  entrance  was  the 
only  zone  that  demonstrated  any  distinct  change  in  plan  to  base  sedimentation. 
These  variations  were  rather  subtle,  less  than  ±10  percent. 

39.  Plans  A,  B,  and  C,  involving  extensions  to  the  north,  resulted  in 
reduced  sedimentation  within  the  Nansemond  zone.  The  largest  change,  about  a 
9  percent  reduction  in  sedimentation,  was  indicated  for  plan  C,  the 
alternative  that  involved  only  northward  extension. 

40.  Increased  sedimentation  within  the  Nansemond  zone  was  indicated  for 
plans  D,  E,  and  F,  which  involve  only  westward  extensions  of  Craney  Island. 

The  largest  increases  were  for  plans  D  and  F  (about  4  and  7  percent,  respec¬ 
tively),  Craney  Island  extensions  which  do  not  extend  to  the  mainland. 

41.  Alternate  plans  A,  D,  and  F  may  impact  water  quality  characteris¬ 
tics  as  a  result  of  a  reduced  circulation  zone  between  the  Craney  Island 
extension  and  the  mainland. 

42.  Formation  of  the  circulation  cell  off  Newport  News  and  Hampton 
Flats  will  continue  under  all  of  the  alternatives  examined.  The  general  loca¬ 
tion  of  the  convergence  zone  was  also  unaffected. 


43.  Potential  three-dimensional  circulation  impacts  to  the  reported 
frontal  phenomenon  cannot  be  quantified  with  the  available  two-dimensional 
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APPENDIX  A:  FINITE  ELEMENT  MODELING* 


1.  The  TABS-2  numerical  models  used  in  this  effort  employ  the  finite 
element  method  to  solve  the  governing  equations.  To  help  those  who  are 
unfamiliar  with  the  method  to  better  understand  this  report,  a  brief  descrip¬ 
tion  of  the  method  is  given  here.  For  a  more  thorough  treatment,  see 
Zienkiewicz  (1971)**  or  Desai  (1979). 

2.  The  finite  element  method  approximates  a  solution  to  equations  by 
dividing  the  area  of  interest  into  smaller  subareas,  which  are  called  ele¬ 
ments.  The  dependent  variables  (e.g.,  water-surface  elevations  and  sediment 
concentrations)  are  approximated  over  each  element  by  continuous  functions 
which  interpolate  in  terms  of  unknown  point  (node)  values  of  the  variables. 

An  error,  defined  as  the  deviation  of  the  approximation  solution  from  the 
correct  solution,  is  minimized.  Then,  when  boundary  conditions  are  imposed,  a 
set  of  solvable  simultaneous  equations  is  created.  The  solution  is  smooth  and 
continuous  over  the  area  of  interest. 

3.  In  one-dimensional  problems,  elements  are  line  segments.  In  two- 
dimensional  problems,  the  elements  are  polygons,  usually  either  triangles  or 
quadrilaterals.  Nodes  are  located  on  the  edges  of  elements  and  occasionally 
inside  the  elements  The  interpolating  functions  may  be  linear  or  higher 
order  polynomials.  Figure  A1  illustrates  a  quadrilateral  element  with  eight 
nodes  and  a  linear  solution  surface. 

4.  Most  water  resource  applications  of  the  finite  element  method  use 
the  Galerkin  method  of  weighted  residuals  to  minimize  error.  In  this  method 
the  residual,  the  total  error  between  the  approximate  and  correct  solutions, 
is  weighted  by  a  function  that  is  identical  with  the  interpolating  function 
and  then  minimized.  Minimization  results  in  a  set  of  simultaneous  equations 
in  terms  of  nodal  values  of  the  dependent  variable  (e.g.,  water-surface  ele¬ 
vations  or  sediment  concentration) .  Time-dependent  problems  can  have  the  time 
portion  solved  by  the  finite  element  methods,  but  it  is  generally  more  effi¬ 
cient  to  express  derivatives  with  respect  to  time  in  finite  difference  form. 

5.  The  finite  difference  method,  sometimes  called  FDM,  is  another  tech¬ 
nique  used  in  numerical  modeling.  The  FDM  solves  mathematical  models  by 

*  Paragraphs  L-4  of  this  appendix  from  Stewart,  Daggett,  and  Athow  (1985). 

**  All  references  cited  in  this  appendix  are  listed  in  the  References  at  the 
end  of  the  main  text. 


approximating  derivatives  with  differences  in  the  value  of  variables  over  fi¬ 
nite  intervals  of  space  and  time.  It  requires  discretization  of  space  and 
rime  into  more  or  less  regular  grids  of  computational  points.  The  finite  dif 
ference  method  obtains  solutions  to  approximate  equations. 

6.  In  summary,  the  finite  element  method,  sometimes  called  FEM,  pro¬ 
vides  a  means  of  obtaining  an  approximate  solution  to  a  system  of  governing 
equations;  the  partial  differential  equations  are  transformed  into  finite  ele 
ment  form  and  then  solved  in  a  global  matrix  system.  The  solution  is  smooth 
over  each  element  and  continuous  over  the  computational  network,  and  the  spa¬ 
tial  integral  of  the  error  is  minimized.  In  comparison,  conventional  finite 
difference  methods  provide  a  means  of  obtaining  a  solution  to  approximate 
equations  of  the  governing  equations;  the  partial  differential  equation  terms 
are  usually  replaced  by  difference  quotients  and  solved  at  discrete  points. 

7.  Simply  stated,  conventional  finite  difference  methods  approximate 
the  equations  to  be  solved,  then  provide  solutions  to  the  approximate  equa¬ 
tions  giving  approximate  answers  for  discrete  points.  Finite  element  methods 
approximate  the  form  of  the  solution  and  then  solve  the  governing  equations 
providing  a  continuous  approximate  solution  over  the  modeled  area  of  interest 
Both  finite  difference  and  finite  element  methods  provide  approximate  solu¬ 
tions  to  the  same  basic  equations  for  conservation  of  mass  and  momentum. 
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APPENDIX  B:  THE  HYDRODYNAMIC  MODEL,  RMA-2V* 

1.  The  generalized  computer  program  RMA-2  solves  the  depth-integrated 
equations  of  fluid  mass  and  momentum  conservation  in  two  horizontal  direc¬ 
tions.  The  form  of  the  solved  equations  is 
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horizontal  flow  velocity  in  the  x-direction 
t  ime 

distance  in  the  x-direction  (longitudinal) 
horizontal  flow  velocity  in  the  z-direction 
distance  in  the  z-direction  (lateral) 
acceleration  due  to  gravity 
water  depth 

elevation  of  the  bottom 

normal  turbulent  exchange  coefficient  in  the  x-direction 
fluid  density 

tangential  turbulent  exchange  coefficient  in  the  x-direction 


*  This  appendix  from  Stewart,  Daggett,  and  Athow  (1985). 
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to  =  angular  rate  of  earth's  rotation 
(j>  =  latitude 

C  =  Chezy  roughness  coefficient 

E,  =  coefficient  relating  wind  speed  to  stress  exerted  on  the  fluid 
V  =  wind  velocity 

¥  =  angle  between  wind  direction  and  x-axis 
e  =  tangential  turbulent  exchange  coefficient  in  the  z-direction 
t  ^  =  normal  turbulent  exchange  coefficient  in  the  z-direction 

2.  The  Chezy  roughness  formulation  of  the  original  code  was  modified  in 
the  input  portion  so  that  Manning's  n  roughness  coefficients  may  be  speci¬ 
fied  from  input  Manning's  n  values  and  initial  water  depth. 

3.  Equations  Bl,  B2,  and  B3  are  solved  by  the  finite  element  method 
using  Galerkin  weighted  residuals.  The  elements  may  be  either  quadrilaterals 
or  triangles  and  may  have  curved  (parabolic)  sides.  The  shape  functions  are 
quadratic  for  flow  and  linear  for  depth.  Integration  in  space  is  performed  by 
Gaussian  integration.  Derivatives  in  time  are  replaced  by  a  nonlinear  finite 
difference  approximation.  Variables  are  assumed  to  vary  over  each  time  inter¬ 
val  in  the  form 

/( t )  =  /( 0)  +  at  +  btC  t  <  t  <  t.  (B4) 

o  —  i 

which  is  differentiated  with  respect  to  time,  and  cast  in  finite  difference 
form.  Letters  a  ,  b  ,  and  c  are  constants.  It  has  been  found  by  experi¬ 
ment  that  the  best  value  for  c  is  1.5  (Norton  and  King  1977).* 

4.  The  solution  is  fully  implicit  and  the  set  of  simultaneous  equations 
is  solved  by  Newton-Raphson  iteration.  The  computer  code  executes  the  solu¬ 
tion  by  means  of  a  front-type  solver  that  assembles  a  portion  of  the  matrix 
and  solves  it  before  assembling  the  next  portion  of  the  matrix.  The  front 
solver's  efficiency  is  largely  independent  of  bandwidth  and  thus  does  not  re¬ 
quire  as  much  care  in  formation  of  the  computational  mesh  as  do  traditional 
solvers. 

5.  The  code  RMA-2V  is  based  on  the  earlier  version  RMA-2  (Norton  and 
King  1977)  but  differs  from  it  in  several  ways.  First,  it  is  formulated  in 

*  All  references  cited  in  this  appendix  are  listed  in  the  References  at  the 
end  of  the  main  text. 
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terms  of  velocity  (v)  instead  of  unit  discharge  (vh) ,  which  improves  some 
aspects  of  the  code’s  behavior;  it  permits  drying  and  wetting  of  areas  within 
the  grid;  and  it  permits  specification  of  turbulent  exchange  coefficients  in 
directions  other  than  along  the  x-  and  z-axis. 
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APPENDIX  C:  THE  SEDIMENT  TRANSPORT  MODEL,  STUDH 


1.  The  generalized  computer  program  STUDH  solves  the  depth-integrated 
convection-dispersion  equation  in  two  horizontal  dimensions  for  a  single  sedi¬ 
ment  constituent.  The  form  of  the  solved  equation  is 
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where 

3 

C  =  concentration,  kg/m 
t  =  time,  sec 

u  =  flow  velocity  in  x-direction,  m/sec 
x  =  primary  flow  direction,  m 
w  =  flow  velocity  in  z-direction,  m/sec 
z  =  direction  perpendicular  to  x  ,  m 

2 

D  =  effective  diffusion  coefficient  in  x-direction,  m  /sec 

X  2 

D^  =  effective  diffusion  coefficient  in  z-direction,  m  /sec 

a  =  coefficient  for  the  source  term,  1/sec 

3 

<*2  =  equilibrium  concentration  portion  of  the  source  term,  kg/m  /sec 
STUDH  is  related  to  the  generalized  computer  program  SEDIMENT  II  (Ariathurai, 
MacArthur,  and  Krone  1977)*  developed  at  the  University  of  California,  Davis, 
under  the  direction  of  R.  B.  Krone.  STUDH  is  the  product  of  joint  efforts  of 
WES  personnel  (under  the  direction  of  W.  A.  Thomas)  and  R.  Ariathurai 
(Resource  Management  Associates) . 

2.  The  source/sink  terms  in  Equation  Cl  are  computed  in  routines  that 
treat  the  interaction  of  the  flow  and  the  bed.  Separate  sections  of  the  code 
handle  computations  for  clay  bed  and  sand  bed  problems.  In  the  tests  de¬ 
scribed  here,  only  clay  beds  were  considered.  Equation  Cl  is  solved  by  the 
finite  element  method  using  Galerkin  weighted  residuals.  Like  RMA-2V,  which 
uses  the  same  general  solution  technique,  elements  are  quadrilateral  and  may 
have  parabolic  sides.  Shape  functions  are  quadratic.  Integration  in  space  is 
Gaussian.  Time-stepping  is  performed  by  a  Crank-Nicholson  approach  with  a 


weighting  factor  (theta)  of  0.66.  The  solution  is  fully  implicit  and  a  front- 

type  solver  is  used  similar  to  that  in  RMA-2V. 

3.  Several  options  are  available  for  computing  bed  shear  stress,  t,  , 

b 

using 


2 

Tb  =  PU* 

where 

p  =  water  density 
=  shear  velocity 

The  Manning  form  of  the  shear  stress  equation  was  used  in  this  study 
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CME[D 
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where 


Krone 


u  =  flow  velocity 
n  =  Manning's  roughness  value 
g  =  acceleration  due  to  gravity 

CME  =  coefficient  of  1  for  SI  units  and  1.486  for  non-SI  units 
D  =  flow  depth 

4.  Deposition  rates  for  clay  beds  were  calculated  with  the  equations  of 

(1962) : 
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where 

Vg  =  fall  velocity  of  a  single  particle 
=  critical  shear  stress  for  deposition 
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C  =  critical  concentration  =  300  mg/£ 

v'  -  V  /c4'3 
k  sc 

5.  Erosion  rates  were  computed  by  a  simplification  of  Partheniades 
(1962)  for  particle-by-particle  erosion.  The  source  term  was  computed  by 
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where 

P  =  erosion  rate  constant 

=  critical  shear  stress  for  particle  erosion 
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